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ABSTRACT 

The growth of supermassive black holes and their host galaxies are thought to be linked, but 
the precise nature of this symbiotic relationship is still poorly understood. Both observations and 
simulations of galaxy formation suggest that the energy input from active galactic nuclei (AGN), as 
the central supermassive black hole accretes material and grows, heats the interstellar material and 
suppresses star formation. In this Letter, we show that most host galaxies of moderate-luminosity 
supermassive black holes in the local universe have intermediate optical colors that imply the host 
galaxies are transitioning from star formation to quiescence, the first time this has been shown to be 
true for all AGN independent of obscuration. The intermediate colors suggest that star formation 
in the host galaxies ceased roughly 100 Myr ago. This result indicates that either the AGN are 
very long-lived, accreting for more than 1 Gyr beyond the end of star formation, or there is a ^100 
Myr time delay between the shutdown of star formation and detectable black hole growth. The first 
explanation is unlikely given current estimates for AGN lifetimes, so low-lumiosity AGN must shut 
down star formation before substantial black hole accretion activity is detected. The scarcity of AGN 
host galaxies in the blue cloud reported here challenges scenarios where significant star formation and 
black hole growth are coeval. Lastly, these observations also strongly support the 'Unified Model' of 
AGN as the host galaxy colors are independent of obscuration towards the central engine. 

Subject headings: galaxies: evolution, galaxies: formation, galaxies: active, galaxies: nuclei, X-rays: 
galaxies 



1. INTRODUCTION 

Several works have noted that obscured AGN appear 
to be prevalent in the 'green valley' on the color- 
magnitude diagram, between actively star-forming 
galaxies in the blue cloud and passively evolving 
galaxies on the r e d sequence (ISchawinski et all l2007bl: 
Salim et alj l2007t iNandra et al l l2007t ISilverman et all 
20081: IConstantin et all 120081: iBundv et all 120081: 
Georgakakis et all I2008T ISchawinski et all 12009( 1 The 



presence of AGN in galaxies with such intermediate 
optical colors has been interpreted as evidence for the 
role of AGN in the suppression of star formation, so 
that AGN host galaxies are transitioning from the blue 
cloud to the red sequence. However, the obscuration 
that affects the central nuclear region might also be 
reddening the galaxy colors systematically. Therefore, 
to understand the connection between star formation 
and black hole growth, we need to study a complete, 
unbiased sample of AGN. 

In particular, do all AGN, including both unobscured 
and highly obscured, Compton-thick AGN, lie in the 
green valle y, as one might hypothesize from the Uni- 
fied Model (|Urrv fc Padovanilfl^lAntonucxilflQQl ? In 
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this Letter, we test the universality of the association of 
AGN with 'green' optical colors using a complete sample 
of obscured and unobscured AGN host galaxies at very 
low redshift and then discuss the implications of the ob- 
served optical color distribution for the role of moderate- 
luminosity AGN in the suppression of star formation. We 
include for the first time not just obscured AGN, but also 
unobscured and highly obscured sources, and can thus 
draw strong conclusions on the host galaxies of AGN - 
not just where on the color-magnitude relation they lie, 
but also where they are absent. 

2. SWIFT BAT SAMPLE SELECTION 

In order to select a complete sample of both unob- 
scured and obscured AGN, we use hard X-ray selec- 
tion, which is nearly unbiased with respect to obscu- 
ration, yet remains highly efficient, unlike IR selection 
techniques. The Burst Alert Telescop e (BAT) onboard 
the Swift satellite (jGehrels et al.ll2003 ) observes photons 
in the 14-195 keV range in which even the most highly 
obscured AGN can be detected. 

We start with the sample of AGN found by the Swift 
BA T in the first 9 mont hs of operation across the entire 
sky (|Tueller et al.| [2008). then limit ourselves to objects 
that overlap with the Sloan Digital Sky Survey ( SPSS; 
lYork et alJ 120001 : lAdelman-McCarthv et all I2008T ) and 
have redshifts of 0.01 < z < 0.07, yielding 21 AGN in the 
luminosity range of Li4_i 95 k c v = 10 42 7 — 10 44 5 erg s _1 
with a median of 10 44 erg s _1 . The AGN include the 
entire range from entirely unobscured to highly obscured, 
and have moderate luminosities, and thus accretion rates, 
typical of the local universe. Quasars exhibit even higher 
luminosities and accretion rates, but are rare in the low 
redshift universe; our sample probes almost all of the 
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Fig. 1. — Two examples of the central point source subtraction for SDSS gr images of a morphologically disturbed galaxy (top row) and 
an elliptical galaxy (bottom row). In panels (a), we show the original image. In panel (b), the central point source is subtracted. In panel 
(c), the galaxy model is subtracted. The residuals show the blue central point source, and in the case of the merger galaxy, tidal tails and 
shells. 



significant black hole growth in the local universe. 

3. HOST GALAXY COLORS AND CENTRAL POINT 
SOURCE SUBTRACTION 

To study the host galaxy optical p roperties, we use 
imaging data from the SPSS PR6 (|York et all 120001 
lAdelman-McCarthv et al.ll2008h . The main challenge in 
understanding the host galaxy colors of AGN is the con- 
tamination of the non-stellar central point source of the 
AGN. Peep X-ray surveys can deliver samples of X-ray 
selected AGN down to moderate X-ray luminosities out 
to considerable redshifts, but the contamination by the 
AGN light can make comparing AG N host galaxies to 
their normal counterparts meaningless (|Simmons fc Urrvl 
120081 ). At high redshifts (z > 0.7), it is difficult to sepa- 
rate host galaxy from central source even with the superb 
angular resolution of the Hubble Space Telescope. At very 
low redshifts however, ground-based optical imaging data 
become comparable, even superior to, the best Hubble 
images of high redshift gala xies. 

We use GALFIT v 2.0.3 ijPeng et al.ll2"002h . modified 
by I Simmons k, Urrvl ((2008) for the optimal removal of 
central point sources, to reliably and robustly separate 
the central AGN light from the extended host galaxy. In 
Figure 1, we show two examples of the results of this 
process. We remove from our sample any sources where 
the point source is brighter than the host galaxy and 



sources where the host half-light radius is less than two 
pixels in order to robustly measure the host galaxy colors 
down to ~0.1 mag. This leaves us with 16 sources from 
the Swift BAT/SPSS sample. This criterion introduces 
a slight bias against the most luminous AGN in the BAT 
sample (Lx > 10 44 erg s _1 ). 

4. RESULTS 

4.1. The Optical Colors of Normal Galaxies and Hard 
X-ray selected AGN Hosts 

In order to compare AGN host galaxies to normal 
galaxies, we construct a volume-limited sample of galax- 
ies from SPSS in the same redshift range with M r < — 19 
mag to select galaxies sampling the entire luminosity 
range of the AGN host galaxies. In Figure 2, we show the 
color-magnitude diagram for the normal galaxy popula- 
tion, featuring the color bimodality of the red sequence, 
blue cloud and the green valley in between. On top of 
this, we plot the host galaxy population of hard X-ray 
selected AGN (large filled circles). Obscured (red) and 
unobscured AGN (blue) clearly reside in the green val- 
ley. No AGN of similar hard X-ray luminosity, regardless 
of obscuration lies in the region typical of star-forming 
galaxies in the blue cloud. 

We use a Kolmogorov-Smirnov (KS) test to assess the 
significance of the difference between the host galaxy col- 
ors of the comparison sample and the AGN host galax- 
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Fig. 2. — The color-magnitude diagram for the host galaxies of X-ray-selected AGN (filled points) and a matched sample of normal 
galaxies (contours, doubling at each level). The g — r color is a good tracer of the age of the stellar population dominating the optical light. 
The galaxy comparison sample from SDSS clearly shows the red sequence and blue cloud. After removing the central point source, the 
X-ray selected AGN appear in the green valley in between. The large points are hard X-ray selected AGN detected by Swift BAT, colored 
according to obscuration (log Nn < 22 -blue; logA^u > 22 - red). The small points are the unobscured AGN detected by ROSAT. On the 
right-hand, we plot g — r color histograms of both the galaxy comparison sample and the AGN sample. 



ies and find that they are different at very high signif- 
icance (3> 3cr). In order to strengthen this statement, 
we compare the color distribution in luminosity bins. In 
all four bins of 1 magnitude each between M r = —23 
to M r = — 19, the g — r color distribution of AGN host 
galaxies and of their normal counterparts are found to 
be significantly different at far more than the 3<r level. 

4.2. Locating the Obscuring Material 

The unobscured AGN show the same green host galaxy 
colors as the obscured AGN, meaning the obscuration 
must be located in the nucleus, not distributed across the 
galaxy. A KS test of the g — r color distribution of the ob- 
scured and unobscured BAT sources yields that they are 
consistent with being drawn from the same parent distri- 
bution at the 2a level. The obscuration seen by the AGN 
must therefore be circumnuclear, as in the Unified Model 
of AGN, such that the gal axy colors are not stron gly af- 
fected by the obscuration. iTreister fc Urr v (2006) show 
that there is evidence that this may not be the case at 
higher redshifts, towards the peak of AGN and galaxy 
merger activity, implying evolution of the obscuration. 



4.3. The Optical Colors of Soft X-ray Selected AGN 

Clearly, unobscured AGN appear to have the same 
green host galaxy colors as obscured AGN - but the sam- 
ple is small. To improve statistics, we obtain a larger 
sample of unobscure d AGN using data fr om the Rontgen 
Satellite (ROSAT; IVoges et all 119991: lAnderson etail 
12001 . which is sensitive in the 0.1 to 2.4 keV energy 
range, and thus to mostly unobscured AGN. Restricting 
ourselves to the same SDSS area and redshift range, we 
obtain a sample of 161 unobscured AGN with soft X-ray 
luminosities of £0.1-2.4 kcV of 10 42 — 10 44 3 erg s^ 1 and 
a median of 10 42 9 erg s _1 . We limit our selection to a 
minimum of 10 42 erg s _1 to avoid contamination with 
X-rays from star formation and X-ray binaries. The en- 
ergy ranges observed by Swift and ROSAT are different; 
however, assuming an unobscured AGN spectrum with 
r = 1.7, the median luminosities of our ROSAT and 
Swift BAT samples are comparable within a factor of 3. 
A KS test affirms that they are consistent with being 
drawn from the same parent distribution, but note that 
the Swift BAT sample is small. 

We again subtract the optical central point sources and 
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Fig. 3. — Time scales for 'green' optical colors of AGN host galax ies. We plot three star formation scenarios probing the limits of possible 
color evolution based on stellar evolution models (jMaraston 2005j). Model 1 is a moderate 10% mass-fraction starburst on top of an old 
8 Gyr stellar population. This burst is somewhat extinguished with E(B-V) = 0.1 and declines with an e-folding timescale of 100 Myr. 
Model 2 represents an extreme scenario of a small 5% mass fraction declining effectively instantaneously, with r = 10 Myr, extinguished 
with the maximum allowable E(B-V)= 0.2. Model 3 has a constant star formation history spanning 10 Gyr with E(B-V) of 0.05; such 
galaxies never leave the blue cloud in most cases. That 95% of AGN host galaxies lie above g — r = 0.45 means that Model 3 is ruled out 
definitely. Models 1 and 2 are consistent with our results, but imply implausibly long AGN lifetimes unless there is a significant time delay 
between the suppression of star formation and the onset of detectable AGN activity. 



add the remaining 161 sources with robust colors to Fig- 
ure 2 (small blue points). Like the unobscured AGN 
detected by Swift BAT, the ROSAT-detected AGN have 
green host galaxy optical colors. The larger numbers al- 
low us to probe the bluest and reddest colors of AGN 
host galaxies, and in particular, to see how many ob- 
jects lie in the blue cloud. We see no AGN in the peak 
of the blue star-forming galaxy population (lower right 
of Figure 2). Due to the range of obscuration repre- 
sented in our sample, this statement must hold for all 
local moderate-luminosity AGN. 

Could we be systematically missing accreting black 
holes at the lowest galaxy mass covered by our sample 
at M r ~ —19? At the bluest colors, this corresponds 
to ~ 10 10 M(T). We convert this to a black hole mass 
(jHaring fc Rixj|2004h and assume that it radiates at effi- 
ciencies as low as a few percent of the Eddington limit. 
Such an AGN would emit a soft X-ray luminosity of 
~ 10 42 erg s _1 and thus should have been detected by 
ROSAT. The lack of blue AGN hosts is therefore signif- 
icant and we conclude that the majority of AGN host 
galaxies in the local universe are in the green valley. 

4.4. A High Merger Fraction for X-ray Selected AGN 
Host Galaxies 

After subtracting the galaxy model with GALFIT, 
the residual image contains both the central AGN point 
source, and any remaining morphological disturbances, 
such as tidal tails, spiral arms and shells (e.g. Figure 1 



top). We perform the same galaxy model subtraction to 
a sample of normal galaxies matched to the AGN sam- 
ple in both g — r color and luminosity, and compare the 
incidence of features indicating a recent or in-progress 
merger in both. We find a statistically highly significant 
excess of such features in the AGN host galaxy sample 
over the control sample (43% for BAT source s, 21% for 
contr ol). We further measure the asymmetry (|Conselicd 
2003) in both the g- and r-band images and find that 
the AGN host galaxies are significantly more asymmet- 
ric than those in the control sample. This observation 
compares favorably to models in which AGN phases are 
driven by mergers (e.g. iHopkins et aT1l2006h . 

5. STELLAR POPULATION ANALYSIS 

We can use stellar population synthesis as a clock to 
interpret the green valley optical colors. In Figure 3, 
we plot g — r colors of three representative star forma- 
tion histories as a fu nction of time, ba sed on the stellar 
population models of iMarastonl (|2005f ). Models 1 and 2 
trace the evolution of realistic and extreme star forma- 
tion shutdown scenarios respectively, and Model 3 follows 
a constant star formation rate. Varying the model pa- 
rameters such as mass-fraction, dust and timescale r fills 
in the parameter space between Models 2 and 3, cluster- 
ing around the 'normal' Model 1. 

A particular concern is whether we can account for the 
green color of AGN host galaxies by invoking substantial 
ongoing, but obscured, star formation. Any extinction 
seen by the host galaxy would also be seen by the cen- 
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Fig. 4. — In this schematic diagram, we show the black hole accretion rate, star formation rate and optical color as a function of time, for 
the case of rapid (t = 100 Myr) suppression of star formation. In the top panel, we show three Scenarios for the black hole accretion rate. 
Scenario A (red, solid) represents a short-lived AGN/quasar phase that suppresses star formation instantly. This kind of scenario would 
predict a large number of AGN in the blue cloud, and so is ruled out by our observations. In Scenario B (green, dotted), the AGN quickly 
switches on, suppresses star formation and remains at high luminosity for ~ 1 Gyr to account for the high fraction of AGN in the green 
valley compared to the blue cloud. In Scenario C (blue, dashed), the moderate luminosity AGN (like those we detect in the X-rays) does 
suppress star formation, but peaks in luminosity several 100 Myr after the suppression, naturally giving ri se to the obse rved distribution 
of AGN host galaxy colors. Such large time delays have been seen in some cosmological simulations (e.g., Li ct al. 2007), although these 
are for high-luminosity quasars at very high redshift. 



tral engine. The unobscured AGN detected with ROSAT 
all have column densities of Ah < 10 22 cm~ 2 . Convert- 
ing this column density to a dust extinction using the 
observed ratio of iBohlin et~aT1 (|1978[ ). A Y = 0.58(^fr), 
yields a maximum value for the optical extinction E(B- 
V) of ~ 0.2 for a typical i?v =3.1. Substantial ongoing 
star formation requires sufficient amounts of dust, but 
such high levels of dust obscuration cannot be present in 
our AGN host galaxies. 

With Model 2, we put a lower bound on the tran- 
sition time from the blue cloud to the green valley. 
Longer suppression timescales directly result in longer 
time delays to the green valley. In Model 1, the r of 
100 Myr is comparable to the dynamical timescale of 
massive galaxies and thus sets a natural timescale that 
compares favorably with observed shutdown timescales 



inferred in post-starbu rst galaxies (e.g., iKavirai et al.l 
I2007al iWild et alj|200l and for m olecular gas reservoir 
destruction (Schawinski et al. 2003). The extreme Model 
2 assumes virtually instantaneous suppresion (t = 10 
Myr). 

The mass fraction of young stars is also a factor in es- 
timating the transition time. The mass-fraction of young 
stars for the extreme Model 2 of 5% is chosen to reflect 
a minor star burst; smaller star burst at the sub-1% level 
are widespread, b ut fail to move red galaxies off the opti- 
cal red sequen ce (lYi et al.H200llSchawinski et alJl2007at 
IKavirai et alj l20"07b| ). While there are still parameter 
choice that might give green optical colors at very young 
ages, such fine tuning seems implausible for our sam- 
ple of almost 200 objects. The mass-fractions chosen for 
the extreme (5%) and realistic (10%) Models are com- 
mensurate to those in green valley early-type galaxies 
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transitioning from the b lue cloud to the red sequence 
([Schawinski et alJl2007bh . 

For the extreme Model 2, the g — r colors of the AGN 
host galaxies restrict the minimum time elapsed since the 
suppression event to ^60 Myr, as only the bluest 5% of 
the population exhibit colors that correspond to younger 
ages. For the more realistic Model 1, the bluest 5% imply 
typical ages of a few hundred Myr. 

6. DISCUSSION 

The fact that we detect very few AGN in the peak of 
the blue cloud seems to rule out scenarios where a short- 
lived, luminous AGN/quasar phase episode suppresses 
star formation on short timescales (i.e. Scenario A (red) 
in Figure 4). 

Two scenarios instead can account for the lack of 
young, blue AGN host galaxies. The first is that the 
AGN suppress star formation in their hosts, but the typ- 
ical AGN lifetime is long compared to the suppression 
timescale, so that most AGN have evolved away from the 
blue cloud (Scenario B (green) in Figure [3}. Roughly 5% 
of our AGN still have ages less than ~60 Myr, in which 
case Scenario B implies a lifetime of > 1.2 Gyr. However, 
the estimates of lMarconi et al.l (|2004l ) reach >1 Gyr only 
for the most massive black holes. The current consensus 
in the literature (e.g. ,lMartini fc Weinbergll200"ll iMartinil 
12004 iMarconi et all 120041 ) give typical lifetimes of only 
a few 10 8 years, which is smaller than the minimum life- 
time implied by the dearth of blue host galaxies. Such 
lifetime estimates are generally given for accretion phases 
with high Ed dington ratios , with the esimate of a few 
10 s years of IMarconi et al.l (|2004h considering efficien- 
cies down to 0.1 -^Edd- The median Swift BAT source in 
our sample has efficiencies of L ~ O.lLEdd with a large 
scatter if we compute L/L-Edd, for black hole masses es- 
timated from the M, - a relation (|Ferrarese fe Merrittl 
l2000t iGebhardt et al.ll200Q[ ) and apply a bolometric cor- 
rection from Lx to Lboi of 20 to 100. 

The other interpretation (Scenario C (blue) in Figure 
4) , is that star formation is suppressed immediately, but 
the AGN luminosity is detectable only with a delay of 
about 100 Myr. If so, low-luminosity AGN during the 
rising phase, below ~ 10 42 erg s _1 , must be sufficient to 
shut down star formation, as suggested by the destruc- 
tion o f cold molecular gas reservoirs by low-luminosit y 
AGN (|Schawinski et al.l l2007bl : ISchawinski et al.ll2009D . 
Then the optical colors of the green valley AGN are a 
reflection of the time needed for substantial amounts of 
material to lose enough angular momentum to reach the 
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black hole and thus increase the AGN luminosity. This 
time gap of ~ 100 Myr cannot be accounted for by a 
heavily obscured AGN phase, as Swift BAT would have 
detected them. In this Scenario C, the low- luminosity 
phase and not the X-ray bright phase is the cause of the 
shutdown of star formation. Perhaps the low-luminosity 
phase, while radiatively inefficient, is channeling substan- 
tial amounts of kinetic energy into the ISM, similar to the 
'radio' or 'maint enance 'mode invoked in current semi- 
analvtic models dBowe r et al.l 120061: 1 Croton et "ail 120061: 
ISomerville et al.ll2008| ). This Scenario also leaves the op- 
tion that another process causes both the shutdown of 
star formation and the AGN phase with a time delay, 
but such a scenario raises more problems than it solves. 

Our observational result challenges recent results 
from sim ulations of AGN fee d back, at least at low 
redshifts dSpringel et all 120051: iDi Matteo et ail 120051 : 
iHopkins et al.l 120061 : Uohansson et al.ll2008l ). which sug- 
gest that the epoch of maximum black hole growth is co- 
eval with the peak of star formation. The coeval scenario 
(i.e. Scenario A) where the AGN luminosity peaks during 
the suppression phase, would predict a large population 
of high-luminosity AGN in the blue cloud, where we ob- 
serve very few. We note however that some of these sim- 
ulations pertain to the high redshift universe, and that 
the physical processes at low redshift may be different. 
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